Experimental observation of bias-dependent non-local Andreev 

reflection 



S. Russo, M. Kroug, T. M. Klapwijk, and A. F. Morpurgo 

Kavli Institute of Nanoscience, Delft University of Technology, 
Lorentzweg 1, 2628 CJ Delft, The Netherlands 
(Dated: February 2, 2008) 

Abstract 

We investigate transport through hybrid structures consisting of two normal metal leads con- 
nected via tunnel barriers to one common superconducting electrode. We find clear evidence for the 
occurrence of non-local Andreev reflection and elastic cotunneling through superconductor when 
the separation of the tunnel barrier is comparable to the superconducting coherence length. The 
probability of the two processes is energy dependent, with elastic cotunneling dominating at low 
energy and non-local Andreev reflection at higher energies. The energy scale of the crossover is 
found to be the Thouless energy of the superconductor, which indicates the phase coherence of the 
processes. Our results are relevant for the realization of recently proposed entangler devices. 

PACS numbers: 74.45. +c, 74.78.Na, 73.23.-b 
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Andreev reflection (AR) is a well-known process that enables charge transfer across an 
interface between a normal metal and a superconductor Q]. At this interface, an incoming 
electron in the normal metal pairs with a second electron to enter the superconductor, re- 
sulting in a reflected hole. Past work has focused on the case of holes that are reflected 
back into the same electrode from which the incoming electrons originate. However, recent 
theoretical studies have considered the possibility that holes are reflected into a second, 
spatially separated electrode 0, S, 0, ^, Q, S Is, Is, IQ, III. It was shown that this "non-local 
AR" process is equivalent to injecting two spin-ent ang led electrons forming the singlet state 
of a Cooper pair into two different normal leads In this way, non-local AR enables 
the realization of solid-state entanglers 1J| -electronic devices capable of sourcing entangled 
pairs of electrons into nano-electronic circuits- that are of interest for quantum information 
processing. 

One way to investigate the occurrence of non-local AR relies on the following idea. Two 
normal metal electrodes are connected via two tunnel barriers (junctions Jl and J2) to one 
common superconducting electrode. If the separation of Jl and J2 is comparable to the 
superconducting coherence length ^, an electron injected at energy E < A from the normal 
electrode of Jl can propagate as an evanescent wave through the superconductor and pair 
with an electron in the normal electrode of J2|3]. This process results in a hole "reflected" 
into the second electrode, i.e. non-local AR. As holes have the opposite charge of electrons, 
holes undergoing non-local AR generate a voltage difference across J2 that has a sign oppo- 
site to that observed when the superconductor is in the normal state (T > T^f). Therefore, 
in principle, the detection of non-local AR is straightforward: Jl is used to inject current 
into the superconductor and J2 is used as a voltage probe to detect a voltage of the correct 
sign. 

In practice, the situation is complicated by the occurrence of a second process competing 
with non-local AR: electrons iniected from Jl can be transmitted into J2 without being con- 
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verted into holes. This process is known as elastic cotunneling (EC) |^ |JJ| and contributes 
to generate a voltage across J2 that has the same sign as that observed when the supercon- 
ductor is in the normal state. Thus, the sign of the voltage measured across J2 depends on 
whether EC or non-local AR occurs with larger probability. The voltage measured at J2 
may also vanish, if cotunneling and non-local AR occur with exactly the same probability 
for all energies of the injected electrons. As some recent calculations predict [4, 11 1 that this 
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could in fact happen, it is not possible to anticipate which signal -if any- will be measured 
experimentally. For instance, in a recent experiment in which two ferromagnetic leads were 
used as normal electrodes, only the sign corresponding to EC has been observed 12 1 . 
In this paper we report a clear experimental evidence for both non-local AR and elastic 
cotunneling using the experimental strategy just outlined. We find that the magnitude and 
the sign of the measured non-local voltage depend on the bias across the injecting junction. 
At low bias, the observed sign is the same as when the superconductor is in the normal 
state, indicating that EC dominates. At higher bias the sign of the voltage is reversed, 
which indicates the occurrence of non-local AR. The energy scale on which the sign-reversal 
takes place corresponds to the Thouless energy of the superconducting layer. From this we 
conclude that the sub-gap microscopic processes of conduction, non-local AR and EC, are 
phase-coherent . 

A schematic representation of the devices used in our experiments is shown in Fig. CJi. 
The structure is implemented in a Nb/Al multilayer sputtered on a thermally oxidized Si 
substrate using conventional Nb/Al technology [l^. The multilayer consists of two normal 
metal layers (Nl and N2, 50nm Al layers) connected via two tunnel barriers to one common 
superconducor (S). Junction Jl is obtained by in-situ oxidation of the Nl layer and subse- 
quent deposition of Nb. Next, a thin (~ 5 nm) Al layer is sputtered on top of the Nb and 
oxidized in-situ. Finally the top Al layer (N2) is deposited to form junction J2. 
The fabrication process used to pattern the multilayer relies on conventional photolithogra- 
phy combined with chlorine-based reactive ion etching. The junctions area is approximately 
4x8 yum^. Independent electrical connections to the three layers are formed by deposition 
of a 200 nm thick Al/Nb layer on a Si02 mask followed by dry etching. We have checked the 
quality of the tunnel junctions by fitting the differential conductance with the usual BCS 
theory and found that the tunnelling characteristics of junctions Jl and J2 do not show 
any substantial difference. This indicates that the superconducting properties of the Nb/Al 
layer (S) are uniform across its thickness. 

In our devices the separation between the two tunnel barriers is determined by the thickness 
of the S layer, which can be controlled on the nanometer scale. This is crucial, since the 
separation of the tunnel barrier has to be comparable to the superconducting coherence 
length in S, ^ ^ V^oh = 10 — 15 nm ,16] (where /g = "^D/vj ~ 2 nm is the elastic 
mean free path, the diffusion constant D = 1.6 cm? / s and = Jivf/ttA). An optical 
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microscope image of one of our devices is shown in FigEJi. 

All the measurements were performed at T = 1.6 K or higher, with the aluminum elec- 
trodes Nl and N2 in the normal state {T^'' ~ 1.2 K). In the experiment we send current 
through one of the junctions (e.g., Jl) and measure the non-local voltage V"''' across the 
other junction (J2), while maintaining the superconductor at ground. The current bias has 
a dc component and an ac modulation with an amplitude of l/iA at 19.3Hz, and a lock-in 
technique is used to measure the ac component of the non-local signal. This corresponds 
to measuring the contribution given to the non-local voltage by only those electrons which 
have an energy E = eVdc, where Vdc is the dc voltage across Jl. 

Fig. 121 shows the V^J measured as a function of Vdc at two different temperatures (above and 
below Tc), on a sample in which the superconducting layer is 15 nm thick (approximately 
equal to ^). At T = 22.5 K, when the Nb is in the normal state, the sample can be simply 
thought of as a resistance network: the measured signal is large, because of the resistance 
of the thin Nb layer, and bias independent. Microscopically, the signal is due to electrons 
injected into the Nb that have a large probability to diffuse into the lead used as a voltage 
probe. At 1.6 K the Nb is superconducting and the Al in the leads is in the normal state. 
Now the non-local voltage is much smaller and it depends on Vdc- Specifically, V^J reverses 
its sign at Vdc = 270/iV and eventually vanishes at Vdc — 700 fiY, thus on a bias range 
much smaller than the superconducting gap (900 /iV, see Fig. Eb)- 

To investigate if this signal originates from evanescent waves propagating below the super- 
conducting gap, we have measured the non-local voltage in samples with different thickness 
d of the superconducting layer. Fig. IHl compares the data measured in three samples where 
d = 15, 50, and 200 nm, respectively. For the 50 nm sample, a non-local signal reversing 
sign with increasing dc bias is still visible at a bias range much smaller than the supercon- 
ducting gap. However, the magnitude of the signal is approximately 20 times smaller than 
for the sample with d = 15 nm. For the sample with a 200 nm thick superconducting 
layer, no non-local signal is observed. These observations indicate that V^J is very rapidly 
suppressed with increasing the thickness of the superconductor, as expected for evanescent 
waves. 

The comparison of different samples additionally shows that the energy scale on which the 
non-local signal reverses its sign (and eventually disappears) becomes smaller for a larger 
separation of the tunnel barriers. For the d = 15 nm sample the zero crossing energy is 
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~ 300 /ieV and for the d = 50 nm it is ~ 50 fieV (see Fig. These values corre- 
spond well to the Thouless energy Et = hD/cP of the superconducting layers, equal to 
Et — 450 /iV and to E^ — 45 /iV for the d = 15 nm and the d = 50 respectively. 
The fact that the Thouless energy determines the behavior of V^^, indicates that the signal 
originates from quantum-mechanically phase coherent processes. This is to be expected, 
since the transit time Ttr of electrons injected from Jl and transmitted into J2 -as electrons 
or holes- is Ttr — d^ / D ~ 1 — 10 ps, much smaller than the inelastic electron-phonon 



{jph — 1 ns at 1 K in Nb) and electron-electron (xge — 0.1 ns) interaction times[l8|. 
Finding that Et is the relevant energy scale in our measurements also gives an indication 
that non-equilibrium effects jl^ in the superconductor do not play a relevant role in deter- 
mining the behavior of V^^. In fact, these effects depend on the quasiparticle injection rate 
and relaxation times, whose energy dependence is not strongly influenced by phase coher- 
ent propagation in the superconductor (and thus by Et). Note also that non-equilibrium 
effect normally become more relevant at higher bias voltage (when the amount of injected 
charge is larger), whereas the amplitude of the signal V^^ is maximum at Vdc = V and 
vanishes for Vdc well below the gap. The absence of non-equilibrium is consistent with the 
low transparency of our tunnel barriers (T 10^^) and with the fact that quasiparticles 
are injected with energies well below the superconducting gap. In contrast to quasiparticles 
occupying states above A, which may have very long relaxation times, quasiparticles with 
E < A decay very rapidly on the scale of h/A. 

Having established the absence of significant non-equilibrium effects, we conclude that the 
measured non-local voltage V^J is due to phase coherent elastic cotunneling and non-local 
AR. EC is predominant at low bias whereas non-local AR dominates at higher bias, where 
the sign of V^J is negative. That the effect is large and present in all samples (approximately 
10 samples with d=15nm and 50nm have been studied) demonstrates that the sign reverse 
in the non-local voltage is not just a sample-specific effect, as has been observed in InAs/Nb 
structures [2^. 

The measured temperature and magnetic field dependence of V^J (see Fig. 4) are consistent 
with this interpretation. V^J increases with lowering T similarly to what one would expect 
from the convolution of a thermally smeared Fermi distribution with an energy dependent 
transmission probability (and excluding the possibility that the signal is due to quasiparticle 
propagating above the gap). The signal is suppressed by a magnetic field applied parallel 
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to the superconducting layer at B ^ 0.5 T, which is much smaller than the critical field of 
our S layer (higher than 6 T 2^). Since A is only slightly reduced (few percent) by such a 
field, we believe that the main effect of B is the breaking of time reversal symmetry for the 
electron-hole wave injected into the superconductor. Note, however, that at 0.5 T the mag- 
netic flux enclosed by typical electron-hole trajectories in the superconductor {d is smaller 
than the magnetic penetration length in Nb for all samples) is only approximately 0.2 x 0o- 
Our observation of a non-local signal shows that the cancellation of the contribution to V^"' 
due to non-local AR and EC does not occur in the samples investigated here. This cancel- 
lation was theoretically found in models that neglect the effect of Coulomb interaction ^1, 
whereas calculations made for different systems in which interactions in the leads play a 
relevant role did all predict the occurrence of visible effects. Since the effect of 

Coulomb interaction on electronic transport is visible in large-area tunnel junctions of size 



comparable to ours 



21[, we believe that Coulomb interaction may also be relevant here. A 



quantitative interpretation of our experimental results will require the analysis of theoreti- 
cal models more sophisticated than those considered until now, which may have to address 
aspects of our samples that have not been considered so far (e.g. a gradient in the phase of 
the superconducting order parameter, or a small sub-gap density of states induced by the 
presence of the normal electrodes). 

In conclusion, we have reported clear experimental evidence for the occurrence of non-local 
Andreev reflection and elastic cotunneling through a superconducting layer. Our results 
show hat these processes are phase coherent and strongly depend on the energy of the in- 
jected electrons. These flndings are relevant for recent theoretical proposals of quantum 
entangler devices that aim at injecting into two spatially separated normal metal leads the 
spin-entangled electrons forming a Cooper pair. In this context, the energy dependence of 
the probability for non-local Andreev reflection may provide a new way to control the output 
of these entanglers. 

The authors acknowledge helpful discussions with Y. V. Nazarov and P. Samuelsson. This 
work was flnancially supported by NWO/FOM and by NOVA. The work of AFM is part of 
the NWO Vernieuwingsimpuls 2000 program. 
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FIG. 1: a) Schematic cross-section of our sample (not to scale). Two normal electrodes (Nl and N2) 
are connected to a superconducting layer (S) via two tunnel barriers (Jl and J2), whose separation 
d is defined by the thickness of the superconducting layer. The concept of the measurement 
configuration is shown in b): current is injected through Jl and the non-local voltage is measured 
across J2. c) illustrates the non-local AR process: an incoming electron from Nl is transmitted 
as a hole into N2 while a Cooper pair condenses in S. d) Optical microscope image of one of our 
samples (top view). The rectangle in the center is where Jl and J2 are located; Nl, N2 and S label 
the electrical contacts to the respective metallic layers. 
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COTUNNBUm NON-LOCAL AR 




FIG. 2: The non-local voltage V^^ measured across J2, on a device with a d = 15 nm thickness 
of the superconducting layer, for two different temperatures. The upper curve is measured at 
T = 22.5 K -well above T^- and shows a bias-independent non-local voltage due to electrons. At 
1.6K (below ), the non-local voltage is much smaller and depends on the bias Vdc across Jl. At 
low bias, V^^ has the same sign measured in the normal state, indicating that elastic cotunnelling 
dominates. At higher bias, the sign of V^"' is reversed, which indicates the occurrence of non-local 
AR. 



9 




-3-2-10123 
\Wmv) 



FIG. 3: (a) Non-local voltage V^J measured at T = 1.6 K on three samples with different thickness 
of the superconducting layer (d = 15, 50, 200 nm, with a normal state resistance of 4.8, 1.7, and 
0.9 Q respectively). Panels b) and c) show the tunnelling characteristics of junctions, measured in 
two devices with d = 15 and 50 nm respectively. The solid line is a fit based on the BCS density 
of states and shows that good agreement is found with A = 0.9 and 1.45 mV for the two different 
thicknesses of the Nb layer 12 ] . 



these thin superconducting films 3| 



he suppression of the gap in the d = 15 nm sample is typical of 
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FIG. 4: (a) Temperature and (b) magnetic field dependence of the non-local voltage V^^ measured 
as a function of V^c, on a sample with d = 15 nm. 
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